Accessing higher order correlations by time-multiplexing 
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We experimentally measured higher order normalized correlation functions (nCF) of pulsed light 
with a time-multiplexing-detector. We demonstrate excellent performance of our device by verifying 
unity valued nCF up to the eighth order for coherent light, and factorial dependence of the nCF for 
pseudothermal light. We applied our measurement technique to a type-II parametric downconversion 
source to investigate mutual two-mode correlation properties and ascertain nonclassicality. 
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Many optical effects involving coherence phenomena 
like interference, diffraction or radiation from fluctuating 
sources can be related to the concept of mutual correla- 
tion Correlation functions unravel their full poten- 
tial in the theoretical framework of quantum optics, pio- 
neered by the seminal work of Glauber [2]. The absence 
of second-order correlations in the famous antibunching 
experiment for light by Kimble et. al. [3] has rigorously 
demonstrated the corpuscular theory of the photon for 
the first time. Normally ordered correlation functions 
constitute a special case of arbitrary operators in quan- 
tum optics, e.g. density matrices, which can be expressed 
by normally ordered boson operator moments Only 
recently a scheme has been proposed (Hf to connect the 
measurement of correlation functions with these gener- 
alized operator moments, rendering a new approach to 
quantum state characterization. 

The experimental study of correlation functions has 
a long standing history. Intensity correlation measure- 
ments were first performed by Hanbury Brown and Twiss 
(if to determine the diameter of stars. Coherent laser 
beams have been studied by Chopra et al. |7|] who exam- 
ined correlation functions up to the third order. For para- 
metric downconversion (PDC) correlations between sig- 
nal and idler modes have been analyzed up to the fourth 
order [H and nonclassical signs in second-order correla- 
tion functions have been demonstrated for heralded sin- 
gle photons [sl , commonly emplo ying multiple avalanche 
photodiodes (APD) as detectors [lO|. Another emerging 
topic is the phenomenon of ghost imaging, tightly related 
to classical versus quantum correlations 12|. 

Though various experiments have been performed 
along the route of measuring normalized correlation func- 
tions (nCF) g^^^ it remained challenging to measure to 
orders higher than four, especially for nonclassical state 
preparation. In this Letter we will present a scalable and 
integrated scheme to determine higher-order nCFs based 
on time multiplexed detection (TMD) [l3| for pulsed 
sources. 

A normalized correlation function (nCF) g^^^ in gen- 
eral is a time-dependent function of the electromag- 
netic field or, in the context of quantum mechanics, 
of the creation and annihilation operators: g^^^ = 



(aT(ti)---aT(t^)a(t^)---a(^i)) 
(at(tOa(ti)>-(at(^^)a(^^))- 



This definition of normalization 



renders ^^^^ for monochromatic states time-independent, 
a finding which remarkably also applies for single-mode 
pulsed ensemble measurements. In fact, any time- 
dependence of g^^^ can be directly related to multimod- 
eness, such that its measurement yields information for 
distinguishing single-mode from multimode states. Con- 
ventionally, the time dependence of nCF is probed in the 
regime of CW light with help of detectors much faster 
than the amplitude fluctuations. In opposition, we focus 
here on the measurement of g^^^ for pulsed sources [8[ 
with a pulse envelope much shorter than the detector re- 
sponse time. An intuitive physical interpretation is that 
the integration over a deltalike intensity pulse essentially 
samples the nCF only at the maximum of the pulse en- 
velope. Any multimodeness can then either be described 
as a temporal substructure within the nCF of the pulse, 
or more conveniently in an effective single- mode picture: 
The photon number statistics, for instance, of a thermal 
light source converges to a Poissonian distribution with 
increasing number of modes, and is caused by a convo- 
lution effect Q. This comes along with a decrease of 
g{n) for the described scenario. Note, however, that this 
decrease is not due to temporal smoothing but in our con- 
text of ultrashort pulses is only caused by the multimode 
character. 

Other features of g^'^^ also provide an easy way for 
classifying quantum states: The second-order nCF is 
a unique operational measure for non-classicality [HI 
if ^/(2)(o) < 1. For various applications like quantum 
key distribution or Schrodinger cat state preparation, 
great interest lies in the feasibility to prepare m-photon 
states [l^. Here, correlation functions can be used 
to quantify the fidelity of an m-photon state since all 
g^^^ — ^^^-'^) must be zero for n > m. Lately, 
a profound and highly applicable theoretical framework 
has been introduced [17] to distinguish classical radiation 
from nonclassical one via correlation properties. For ex- 
perimentalists an intriguing feature of a nCF is its loss 
independence^ i.e no a priori assumptions about channel 
or detector efficiency are required. 

The theoretical model of our experimental setup corre- 
sponds to a hierarchical beamsplitter cascade of stages 
as shown in Fig.[TJ Losses can be accounted for by adding 
a virtual beamsplitter in front of each detector, such that 
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only the transmitted amplitude is sensed. In the follow- 
ing we mathematically elucidate this experimental config- 
uration to be harnessed for measuring g'^^^ up to n = 2^. 
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FIG. 1. Spatial mode multiport: Beamsplitter cascade mixing 
signal input with vacuum at each beamsplitter 

Suppose that each beam splitter at stage k transforms 
the input mode sl^-ij and vacuum mode h^-ij (not 



shown) into a/c^2j-i and a/c^2j as illustrated in Fig.[Tl We 
introduce the parameter Sk^i in the transformation for a 
beamsplitter, denoting a transmission amplitude for the 
second index / being odd, and reflection amplitude for / 
being even: 



a/c,2j = Sk,2j^k-l,j 



- ^ (1) 



Starting from mode ao,i the signal transverses the cas- 
cade to a marked detector j along a unique path 
T^j ^ T^j ^ ^ T^j^' Each mode along this 
route is denoted as kpk. The intensity correlation 

(ajy ^•^aAr,ji • • • ajy ^-^aATjv,) between n distinct and in- 
dependent detector modes can be reordered to read 
• • • di]yj^diNj^ • • • aATji). This expectation value 
can be stepwise expressed by modes of the previous 
beamsplitter stage with help of Eq. [TJ By explicitly trac- 
ing out the vacuum modes and, iterating back to the 
origin, we obtain: 
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The value of g^^"^ for the signal of interest can therefore be 
obtained when dividing the mutual intensity correlation 
of n detectors divided by the intensity of each individual 
detector. All losses, transmission and reflection ampli- 
tudes resembled by the product of Si^k factorize in both 
numerator and denominator and thus cancel out. Hence, 
our calculations confirm ^^^^ to be loss independent. 

(^0,1^0,1)^ (aj^^^.^aATj,) • • • (aj^^^.^aAT,^-^) 

Note that this result can easily be generalized to cross- 
correlate ^(^'^) for two spatial modes a and b with in- 
dependent cascades by considering coincidences of m de- 
tectors in mode a with n detectors in mode b, followed 
by an analogous normalization. 

The measurement of g^^^ yet requires linear detectors 
with respect to the intensity n = a^a. Using linear pho- 
todiodes would pose a straightforward approach. How- 
ever, their signal at the single photon level is concealed 
by the readout and amplification noise. Conventionally, 
APDs are employed that operate in Geiger mode but only 
yield a binary click response upon the detection of pho- 
tons. Since a detection process causes an avalanche in 
the detector's band gap, the number of impinging pho- 
tons cannot be resolved. On the other hand, the proba- 



bility for a click is linear under certain conditions. The 
probability to excite m photoelectrons is given by the 
operator = : ^(r^n)^ exp(— r^n) :, with detection ef- 
ficiency T] and : : denoting normal ordering. Since our 
detector cannot discriminate the number of photoelec- 
trons, any number m greater than one P = Xl^^i Pm 
gives rise to a click. For the regime of states with low click 
probability, P ?7n is dominated by the linear term. In 
the limit of low click probability an APD acts as a lin- 
ear but lossy detector which can be employed to measure 
the loss invariant nCF in Eq. [2j Neglecting dark counts, 
the approximation to the linear behavior of P becomes 
even better for higher losses. The accuracy of the corre- 
lation measurement then increases, however at the cost 
of statistical precision. 

The implementation of the free space spatial mode 
multiport in Fig. [1] even for a modest amount of stages 
becomes very resource consuming. The same network, 
however, has previously been realized as a stable fiber in- 
tegrated TMD [13], depicted in Fig. [2l Here, the spatial 
modes of the cascade are mapped to modes in the time 
domain. The previous theoretical considerations there- 
fore remain identical. One additional fiber loop between 
two subsequent beamsplitters causes a delay between the 
temporal modes, which has been chosen to be twice the 
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FIG. 2. Time mode multiport: Experiment for correlation 
measurements of coherent and chaotic states by a TMD. For 
abbreviations see text. 
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TABLE I. Ah 28 possible second-order coincidence subsets 
for eight total TMD modes. Displayed are the measurement 
results for g^"^"^ after the first 30 s of measurement time for the 
coherent light scenario, yielding a mean g^'^'^ — 1.002 di 0.007. 



APD dead time. 

Coherent and thermal states provide an ideal testbed 
for nCFs. Firstly, they are easy to prepare, and secondly, 
they have a simple analytic expression: g^^^ = 1 for co- 
herent and g^^^ = n\ for thermal states. Coherent states 
are readily provided by a laser source whereas pseudo 
thermal sources can be mimicked by chaotic states, eas- 
ily obtained from a moving speckle pattern [Isl . We em- 
ployed the coherent state and the chaotic state in order 
to validate the TMD's potential for measuring nCFs. 

Our first experiment was conducted to demonstrate the 
capability of measuring g'^^^ with a TMD, and the exper- 
imental configuration is depicted in Fig. [2j A diode laser 
generated 50 ps pulses of single- mode coherent states at 
810 nm with a rate of 20 kHz. The pulse was heavily at- 
tenuated by a neutral density (ND) filter to the single 
photon level. The photons were coupled into one input 
of a multimode fiber based TMD. This yielded an average 
count rate of 3.3 kHz per detector mode. The electronic 
detector signal was processed by a time-to-digital con- 
verter (TDC). The TDC numbered consecutively all laser 
pulses and recorded the time information of an APD-click 
with a precision of 81 ps relative to the laser trigger, much 
better than the integration time of the APD (Ins). The 
chaotic light was generated by focussing the laser onto a 
rotating ground glass (RGG). This way a moving speckle 
pattern was formed and only an areal portion of less than 
an average speckle was coupled into a single-mode fiber 
based TMD with help of two additional apertures (AP). 
Here, the laser repetition rate was increased to 500 kHz 
and an average count rate of 6.6 kHz observed per TMD 
mode. 

Extracting g^^"^ from the recorded data was performed 
by selecting n arbitrary TMD modes. Let all those possi- 
ble ( ^ j selections of n modes be contained in set C. For 
the particular case of n = 2 the set C consists of all pairs, 
exemplarily shown in Table |T] for a ^^^^ measurement of 
a coherent input state. Each combination c G C provides 
an approximation to g'^^^ when dividing the n-fold co- 
incidence probability by the product of individual single 
click probabilities according to Eq. [2l This gives a string 
of values where the mean reaches a value closest to the 
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TABLE n. Measurements of g^^-* for coherent and chaotic 
light sources. Displayed are the theoretical prediction (Th.) 
in the 1^*, the experimental value (Exp.) in the 2^^, and the 
error estimation (Err.) in the 3^^ row. 



real g^'^\ Some benefits are offered when a coincidence 
measurement does not select all 2^ modes and C contains 
multiple elements. Firstly, part of the signal is diverted 
to other modes not being part of the particular selection. 
Hence, measurement accuracy is increased due to higher 
losses considering c alone. Secondly, the signal is not lost 
within the entire multiport network but is transferred to 
other combinations in C. Therefore, good measurement 
statistics can be maintained by averaging over the results 
of C. Additionally, a rough error estimation can be given 
by the standard deviation. 

Our results for higher order correlations are presented 
in Table [III We find an excellent agreement between 
theoretical prediction and experimental findings for both 
coherent and chaotic states. We have proven with confi- 
dence that the TMD can be utilized for accurate measure- 
ment of nCFs and provides the basis for more advanced 
quantum state characterization. Interestingly we find an 
overestimation of the error for n > 5. Contribution to the 
growing standard deviation originates from an increasing 
number of possible combination in C. Another, more fun- 
damental effect is the growth in the standard deviation 
of g'^^^ for chaotic light and large n as a result of strong 
intensity fluctuations. It is this growth that reduces the 
signal-to-noise ratio of multiphoton ghost imaging with 
chaotic light sources [Toj . 

In our second experiment, depicted in Fig. [3^1, we in- 
vestigated correlations and non-classicality between sig- 



4 



nal and idler beams from the nonlinear process of para- 
metric downconversion (PDC). Since the photon number 
between signal and idler strictly match for a PDC state 
1^) — Xln V^l^'^)' measurement of nCF is partic- 
ularly interesting. The beam from an ultrafast (100 fs) 
Ti:Sa laser system, pulsed at 1 MHz at 796 nm) was fre- 
quency doubled in a nonlinear /3-BaB204 (BBO) crystal 
and used to pump a PDC process (order of 100 fs) in 
a periodically poled KTiOP04 (KTP) waveguide chip. 
Any residual pump light was eliminated by spectral fil- 
ters (SF). Signal and idler modes emerged from the type- 
II PDC process with orthogonal polarization and were 
spatially separated by a polarizing beamsplitter. Both 
beams were launched into input fibers of our TMD, 
one being delayed by 400 ns to emulate two independent 
TMDs. 




mean marginal photon number (n) of the PDC state |^) 



lent agreement between measurement and theory without 
any fit parameter, confirming the potential of the TMD 
as well as the Poisson statistics obtained from multimode 
PDC. 

To conclude, we have analyzed in theory and experi- 
ment the use of a beamsplitter cascade for measuring nor- 
malized correlations functions g^^^ with standard APD 
detection. We have implemented the cascade as an in- 
tegrated and easily scalable TMD fiber network. The 
feasibility to correctly measure high order ^^^^ has first 
been ascertained by examining coherent and pseudother- 
mal states of light. After this, we applied our measure- 
ment method to PDC states to study two-mode corre- 
lation functions, the non-classicality of PDC states, as 
well as g^'^^ for heralded single photon states with rem- 
nant higher photon number contributions. We consider 
our experiments to be a first step paving a new way for 
state characterization with TMDs utilizing arbitrary mo- 
ments , which is capable of detecting quantum features 
of correlations independent of detrimental loss effects. 

We would like to thank I. N. Agafonov and W. Vogel 
for fruitful discussions. This work was supported by the 
EU under QAP funded by the 1ST directorate as Con- 
tract no. 015848. 



FIG. 3. a) Experimental setup for mutual correlation mea- 
surements, b) shows the dependency of mutual correlations 
theory (solid) and measurement (squares) for five 
different (fi). c) displays g^'^^ of the idler beam constrained 
on an APD click in the signal beam, hence demonstrating the 
fidelity of heralding a single photon state. 

We parametrized the mutual correlation function 
g,(n,m) _ {b. h h ) ^crms of the mean photon num- 

^ (ata)^(b^b)"^ 

ber for comparing measurement and theory. The model 
assumed our PDC state to exhibit a Poissonian distribu- 
tion pn due to the excitation of many spectral modes [l^ . 
The experiment was conducted with five different sets of 
pump power resulting in different mean marginal photon 
numbers (n). In order to calibrate (n) we applied a loss 
estimation procedure for TMDs as presented in [2o| . The 
results for g^^''^\ g^^'^^ .^^''^^ and ^^^'^^ are presented in 
Fig. [3]3. According to [17| a possible criteria for the vio- 
lation of classicality is g^^''^^ = 7\/ ^(^'^)^(^'^) for 7 > 1. 
Our experimental data clearly affirms this violation with 
a value of 7 ranging between 1.19 and 1.60, depending 
on pump power. We therefore can testify non-classicality 
of 1^), even for lossy detection. Heralding a single pho- 
ton state upon a click in one mode can leave remnants 
of higher order photon contributions in the other mode. 
This can be seen in Fig.[3fc for a conditioned ^^^^ measure- 
ment, proving nCFs to be an extremely powerful tool for 
analyzing state preparation techniques. Note the excel- 
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